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Abstract 

Two dimensional flows on fixed smooth surfaces have been studied in the point of view of 
vorticity dynamics. Firstly, the related deformation theory including kinematics and kinetics 
is developed. Secondly, some primary relations in vorticity dynamics have been extended 
to two dimensional flows on fixed smooth surface through which a theoretical framework 
of vorticity dynamics have been set up, mainly including governing equation of vorticity, 
Lagrange theorem on vorticity, Caswell formula on strain tensor and stream function & 
vorticity algorithm with pressure Possion equation for incompressible flows. The newly 
developed theory is characterized by the appearances of surface curvatures in some primary 
relations and governing equations. 
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V 

£3 ' 1 Introduction 

t/3 | Two dimensional flows/deformations are naturally referred to flows/deformations of con- 

tinuous mediums whose geometrical configurations can be taken as two dimensional surfaces 
embedded in three dimensional Euclidian space. Physically, the values of the thickness of con- 
tinuous mediums are quite smaller than the ones of character scales along flow directions. Such 

Q-f flows/deformations may be divided into two groups. The first one is termed as flows on fixed sur- 

faces that can be considered as models of long wandering rivers on the earth, real floods spreading 
over plains, depressions and valleys, atmosphere motions on planets, creeping flows on curved 
surfaces, flat soap films and so on. The second one is termed as self-motions / deformations as 
models of bubble, cell and capsule deformations, vibrations of solid membranes, oil contamina- 
tions on sea surfaces et al. 

The ways of attaining governing equations of natural conservation laws can also be separated 
into two kinds. The first ki nd is to do app r oxima tions based on general full dimensional govern- 
ing equations. For example, iRoberts &: Lil ( 20061 ) attained models of thin fluid flows with inertia 



on curved substrates based on central manifold approximation. The second kind is to do analysis 
in the point of view of general continuum mechanics. The equations governing fluid motion in 



a surface or interface have been studied by lArisI (J1962I ). Very recently, a novel kind of finite de- 
formation theory of continuous mediums whose geometrica l configurations are two dimensional 
smooth surfaces have been developed bv IXie et al.1 ( 20131 ). It mainly includes the definitions 



of initial and current configurations, deformation gradient tensor with its primary properties, 
deformation descriptions, transport theories and governing equations corresponding to conser- 
vation laws. The general theory is suitable for any kind of two dimensional flows/deformations. 
Particularly, the present s tudy puts focus on flows on fixed smooth surface in the point of view 



of vorticity dynamics (see IWu et al.l . 120051 ) . 



Current Parametric Configuration 




Initial Parametric Configuration 



Figure 1: Sketch of the initial/current and physical/parametric configurations, where n^ 
and Ns denote the surface normal vectors associated to the current and initial physical con- 
figurations respectively, {gi}f =1 and {G^}^ =1 are local co-variant bases. The tangent space 
denoted by TS is span by the co- or contra- variant basis. 



2 Finite Deformation Theory 

2.1 Configurations, Coordinates, Curvatures & Field Operators 

A two dimensional flow on a fixed smooth surface is sketched in Figure [TJ The surface 
is represented by the vector valued mapping S = S(x) e K 3 , and {x l }f =1 is the paramet- 
ric coordinates. The flow/deformation can be represented by the smooth diffeomorphism x = 

o t 

x(£,t) £ ^^(E^, T, x ) in the parametric space, where {x l }f =1 and {£ }\ = \ denote Eulerian and 



Lagrangian coordinates respectively, and £j and S x are termed as the initial and current para- 

o o t t 

metric configurations respectively. Correspondingly, X := S(Sf) and X := S(S :r ) are termed 
as the initial and current physical configurations respectively. The smooth surface can be taken 
as a two dimensional Riemannian manifold with the fundamental quantity of the first kind 
{dij — (9ii9j)}ij=n an d the one of the second kind {bij = (-g^i(x),n^,)}f - =1 . Gaussian curva- 
ture and mean curvature are defined as Kg = det[6y]/det[<7y] and H = b s s respectively. Einstein 
summation convention is adopted throughout the paper with subscripts and superscripts repre- 
senting co- and contra-variant quantities respectively, and quantities with capital/small scripts 
associate with initial/current configurations. 

1 Ft 9 

Generally, the surface gradient operator V = g -Sj is defined as, say <1> G 2? (TS), 
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where o— can be any available algebra tensor operator, Vz denotes the co-variant deriva- 
tive/differentation of the tensor component, e.g. 
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where Fj denotes Christoffel symbol of the second kind. The contra-variant derivative relates 
generally to the co- variant one through V = c i V(. On the other hand, one defines the so 
termed Levi-Civita connection operator V = g l V a 

Dx' 



V o -* = (flf'VjL) o -($?. ff - ® ^) 4 g< o -V a ($?,ff 
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dx' 1 



dx'' 
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V&Jg 



where V_a_ denotes Levi-Civita connection defined on the surface/Riemannian manifold (see 

Qj -' — I I -\ 

Durovin et all 11992 ). It should be noted that the change of the order of co- and contra- variant 

derivatives must be related to Riemannian- Christoffel tensor that can be represented by Gaus- 
sian curvature and metric tensor for two dimensional Riemannian manifolds as revealed by the 
relation 



V p V 9 $?j 



V'Vp^ + K G (6;Si - g sp g iq )^ + K G (g rp g^ - 5 s p 5])& s 



where i? J sp 9 = blfil — b sp b iq denotes the component of Riemannian-Christoffel tensor. 

It is worthy of note that Levi-Civita connection operator is just valid/effective for surface 
tensor fields that just have components on the tangent plane, generally denoted by 2? r (TE) 
where r G N denotes the order. However, the surface gradient operator can be applied to 
arbitrary tensor fields defined on the surface due to its definition is based on differential calculus 
of tensor normed space. 



2.2 Kinematics 

Generally, the velocity of a fluid partial is determined as follows 
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Subsequently, the acceleration is 
A 8V,. , fdV 1 



a - — (f,t) 



(x, t) + V s V s V l ) g t + {b ij ViV 3 ) n =: a l g t + a n n 



dt vs ' ' \ dt 
The material derivative of any surface tensor field is represented by 



(2) 



(3) 



based on the velocity representation ([JJ). The deformation gradient tensor denoted by F is 
introduced by the relation 



£(*(£ + St t),t) - E(x(e, t), t) = F • £(£ + 50 - £(£) , F : 



dx 1 
djA 
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All kinds of deformation descript ions can be repr esented by the deformation gradient tensor 
with its properties as indicated in lXie et al.1 ( 20131 ). 

Similar to the familiar Helmholtz velocity decomposition, one has the relation 



V(x + 5x, t) - V(x, t) = (V ® V) • 5E 



;(VjV i + V i V j )g i ®g j 
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Accompanying (j4|) with the dual relations 

(VjVi - ViVj) =: e ji3 u 3 & 



UJ 



e 3ij V l V j 



(5) 



c 3kl 



the vorticity is defined as u = e VfcV^ns = w ns, where ej^ := [g^g i ,riY\ and e 
[ n E,9 k ,9 l ] are termed as co- and contra- variant components of Eddington tensor that as the 
same to metric tensor are insensitive to co- and contra-variant derivatives. And one has the 
novel velocity decomposition. 

Lemma 1 (Velocity Decomposition) 

V(x + 5x, t) - V(x, t) = D • 5S + ~u x <J£ + V ■ K • «JS, 

where D := |(ViVj + VjVi)g l ® gf- 7 can 6e defined as strain tensor, K := fr^gf* (8> <7 J «s termed 
as surface curvature tensor. The last appended term is due to the contribution of the surface 
curvature. 



2.3 Kinetics 

Basically, so termed intrinsic generalized Stokes formula of the second kind is introduced 
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(6) 



where n := r x n^ is the normal vector of the boundary on the tangent plane in which r is the 
unit tangent vector of the boundary, $ can be any kind of tensor field defined on the surface. It 
is worthy of note that most of the governin g equations for fl ows on surfaces and the ones for thin 
enough shells and plates, as indicated by lArisI (|1962l ) and ISvnge fc Chienl ( 194ll ) respectively, 
ca n be direct l y and readily deduced through ([6|). Its detailed proof can be referred to the paper 
bv lXieet all (J2013h . 

On the mass conservation, it can be represented and analysed 



= J t ^(x,t)da+Ln-(pV)dl 



dp 



hdt 



x, t) do- + l t V 
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where p denotes the surface density. Subsequently, the differential form for mass conservation 
is deduced 



^(x,t)+V-(pV) 
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(7) 



On the momentum conservation, it is considered that the rate of change of the momentum 
are contributed by the surface tension, inner pressure, inner friction and surface force, in other 
words it is keeping valid 
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where F 



ten 



§ t 7 n dl, F™i := -§ t pn dl, F vis := § t pn • [(V^ + V^-) g l ® g j ] dl 
9E 9E 9E 

are represented originally as curve integrals that can be transformed directly to surface inte- 
grals through the intrinsic generalized Stokes formula of the second kind, 7 and p denote the 
coefficients of surface tension and inner friction/viscousity respectively. F sur := J t f sur da is 

represented directly as surface integral, f sur can be the surface densities of fricti o n, gr avity, 
electromagnetic for c e and stochastic force et al. As reviewed by iBoffetta fc Eckel ( 2012 ) and 
Bouchet fc Venailld ( 20121 ) . the surface force usually plays the important role in dynamics of 



two dimensional flows. Furthermore, the component equations of momentum conservation can 
be deduced 



( 



pai = -jj^iP, t) + H [f'ViVjVt + Vi (V S V S ) + K G Vt] + f sur j 
pa n = H{n -p) + fi [2b ij ViVj] + f sur ,n 

where a\ and a n are the components of acceleration on the tangent plane and in the normal 
direction respectively, as shown in ((2J). It is revealed by (jHJ) that Gaussian curvature accompa- 
nying with the tangent velocity takes part in the momentum conservation on the tangent space, 
however mean curvature accompanying with the surface tension and inner pressure does the 
contribution to the momentum conservation in the surface normal direction. 
On the moment of momentum conservation, one has the general conclusion 

Lemma 2 (Moment of momentum Conservation) Generally, the mechanical action im- 
posed on the boundary can be represented by the so termed surface stress tensor 

t = tj9i <8 9 j + t%g { &> n s 

then the moment of momentum conservation takes the form 

g'-(txg,) + m s = 06 M 3 
where mj denotes the surface force couple. 

Proof: Based on the surface stress tensor, the momentum conservation can be represented 
generally as 

A p a da = f t n • t dl + A f s da =^ p a = V • t + / s (9) 

is Jas Jt, 

On the other hand, the moment of momentum conservation is represented as 

A p a x £ dcr = f t (n • t) x £ dl + A / E x E da + A rn^ da 

Jt, JdT, Jt, Jt 



with the differential form 



T 

paxS = V'(txS) + / s xS + m s 



V-t x£ + s'-(txs. 



+ / s x S + m s (10) 



Substituting ([9]) into (fTUj) . the proof is completed. Furthermore, thanks to the representation 

g l • (t X 0,) = -t ij e ij3 riT + Vgi-t 2 ^ 1 + t\g 2 ), g := det[ 5ij ] 
it can be concluded that 

Corollary 1 (On Representations of Surface Stress Tenor) The symmetry of the com- 
ponents of surface stress tensor on the tangent space, i.e. tij = tji, corresponds to the vanishing 
of the component of surface force couple in the surface normal direction. And the appearance of 
surface stress tensor in the surface normal direction, i.e. t l 3 7^ 0, corresponds to the existence 
of components of surface force couple on the tangent space. 

As an application, the stress tensor corresponding to the actions of surface tension, inner 
pressure, inner friction as discussed previously takes the form 

t = ( 1 -p)I + p(V J V i + V i V J )g i ®g j , I := m tf % g j (11) 

that corresponds to the case of full zero surface force couple as usually satisfied by Newtonian 
fluid flows. 



On energy conservation, one has the identity 






j(x,t)da + ^n-( P V)(e+^-j dl 
L(n-t)-Vdl + I / E • V da + I g s da 



where e denotes the surface density of internal energy and q^ the heat flux. Furthermore, using 
@, one can arrive at 
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where the last two identities are due to the adoption of the constitutive relation ([TT 



3 Vorticity Dynamics 

3.1 General Theories 

Firstly, the following identity has been derived 
Lemma 3 



V x (b-F) 



AT E = |F| (V x 6) • n s , VfoeTS, |F| := 
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12) 



where V = G V_s_ and V = g v _a_ are Levi-Civita connection operators, N^ andnj] are sur- 

d£ L dx l 

face normal vectors corresponding to the initial and current physical configurations respectively, 
|F| denotes the determinant ofF, yG:= [G%, G2, JVs]- 



Proof: 
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= |F|(V x b) -n E 
As an application, the governing equation of vorticity can be derived 
Corollary 2 (Vorticity Equation) 

dj 3 = -Ou 3 + (V x a) • n s 
Proof: Let b in the relation (|12p be the velocity V, it reads 



a; 3 = (V x V) • n E 



Vx(V-F) 



iVs 



(13) 



Furthermore, one can do the following deduction 

1 



u) ' 



Vx(y.F) 



N s + 



-0(V x V) ■ n s + (V x a) • n E + V x V 



V x (a • F) + V x ( V • (V 

V\ 2 



s 
V)-F 



AT S 



ns = —Ouj + (V x a) • ns 



where the identities F = (V ® V) • F and d|F|/di = 0|F| fsee lXieet aD . l2013h are utilized. 
As an appendant of the above proof, one has the conclusion 

Corollary 3 (Lagrange Theorem on Vorticity) In the case of the acceleration field is ir- 
rotational, a patch of continuous medium that is initially irrotational will keep irrotational at 
any time, conversely a patch of continuous medium with initially nonzero vorticity will possess 
vorticity at any time although the value can be changed. 

Secondly, the following identity can be readily set up 
Lemma 4 

V x (V x b) = V(V -b) + K G b- Ab, V&eTS, Ab = V • (V x b) 
Proof: 
V x (V x b) = V x \(g p V^) x (btg* 



O.rl' 



giV^_ x [(V p bi)eP i3 n] 
= e 3k "e 3p ^ q (V^)g k = (5 k p 5« - <^)V g (W) fffc = V^V^)^ - V p (W) 5i 
= V t {V k b l )g k - Ab 

Furthermore, one has 

V r {V k V)g k = [V fc (V^) + Ktb s ]g k = V(V • b) + K G {5\5 k - g sl g tk )b s g k 
= V(V -b) + K G b 

where the change of the order of covariant or contra-variant derivatives should be related to 
Riemannian-Christoffel tensor. It is the end of the proof. It should be pointed out that this 
kind of identities is still keeping valid for any tensor field. 

As an application, the governing equation of momentum conservation on the tangent plane 
(18]) can be rewritten as 



>cng l = Vn-/iVxw + 2fiK G V + f sur ,ig l , U := -p + 2/x0 



(14) 



where V x u = e fc3/ ^x(x, t)g L . Subsequently, the following coupling relations can be attained 



dx k 

just by doing the dot and cross products by e on both sides of (fl~4"|) respectively. 

Corollary 4 (Coupling Relations between Directional Derivatives of II and u> ) 

<9IT 

— = pa-e + n(V x w) • e - 2\iK G V • e - / s • e 

oe 



/'■ 



dco 3 
de 



-[pa, e, n s ] + [VII, e, n s ] + [2pK G V, e, n s ] + [/ s , e, n s ] 



for all e 6 TX! s.t |e| = 1. The coupling relations are valid at any point in the flow field. 
Thirdly, the intrinsic decomposition is still valid for any surface tensor, i.e. there exists 



Lemma 5 (Intrinsic Decomposition) 

{e<g>(e- *) - [e, [e,*]] 
, V* G ^ r (TS), eeTS s.i. |e| = 1 
(*-e)<g>e- [[*,e], e] 

where [e, 3>] denotes the cross product e x <1> and so on. 

Proof: Let 3> = 3>?,<7j <8> fl^ S ^ 2 (TS) without lost of generality, and to calculate 

[[*,e], e] = [[& j g i ®g j ,e k g k ], e] = [^e fc e^ 3 ^ ^ n s , e] = $ i j (e fc e 9 )(e^ 3 e, J3 )g i ® g' 
= ^-(e^Xo^f - <^) 5i ® flf' = ($^) ® (e fe g fe ) - $ 

Evidently, one arrives at the conclusion. The other identity can be proved similarly. 

Applied the intrinsic decomposition to the velocity gradient tensor, the strain tensor has the 
following representation, termed still as Caswell formula. 

Corollary 5 (Caswell Formula) On any fixed solid boundary where the general viscous bound- 
ary condition is satisfied by the fluid, namely V = on the boundary, the strain tensor can be 
represented as following 

D = #n (g> n H — (w x n) ® n H — n <g) (u> X n) 

Proof: Basically, a kind of local curvilinear coordinates denoted by {x M ,x^} corresponding to 
the fixed solid boundary can be set up such that the local co- variant basis {g u ,g±} is orthogonal 
on the boundary with g u tangents to the boundary and g± = g^ = n. The details can be 
referred to the appendix. 

Subsequently, the intrinsic decomposition is utilized 

V(g)V = (V<g)V,n)<g)n-[[V<g)V,n],n] 

For the first term on the right hand, one has 

(V (g> V, n) = (V <g) V - V <g> V, n) + (V <g> V, n) = u x n + (V ® V, n) 

where 

(V® V,n) = (g"<8) V _a_V + g- 1 - ® V_a_V,n) = 6n 

V dx" dx 1 - ' 

thanks to the relations 

V^V = Vyg„ + V ll v ± g ± =0gTS 

dx 11 

due to the viscous boundary condition and 

V^_V,n) = (V ± v t, g u + V ± v ± g ± ,g ± ) =V ± v ± = 9-Vy = 6 

For the second term on the right hand, it is naught thanks to 

[V0V,n] = [(V i y)® g ll + (V i y)8 ff 1 , n ] = [(V i y)®g i ,g i ] = (V i y)®[ 9 1 ,g 1 ] = O 

As a summary, one arrives at the relation 

V(&'V = (uxn)(g)n + 6n,(g)n 

Accompanying with its conjugate relation, the proof is completed. 
Readily, one has the following assertion 

Corollary 6 For any two dimensional incompressible viscous flow on any fixed smooth surface, 
on any fixed solid boundary, the directions corresponding to the maximum or minimum rate of 
change of element material arc length with the same absolute value |w 3 |/2 is 7r/4 or Stt/A with 
respect to the tangent direction of the boundary. 



3.2 Particular Theories For Incompressible Flows 

As generally, t he incompressibility for two dimensional flows on fixed smooth surfaces is still 
defined as 6 = 0. IChomaz Sz Cathalaul ( 19901 ) revealed some velocity domains of soap films in 
which they can be considered as two dimensional incompressible flows. Based on the continuity 
equation (J7J), one has 



v,y s 



dx s 
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y/gdx 1 



(^V s )(x,t) = 



,st3 dip 



then the stream function can be introduced through V s = e s j£(x, t). Subsequently, the stream 
function & vorticity algorithm can be derived. 

Lemma 6 (Stream Function &: Vorticity Algorithm for Incompressible Flows) 



AV = 3* 
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V s V s u 3 + 2e^V k (K G V l ) 



+ -e Vfc/surJ 



Proof: The Stream function Possion equation is just from the definition cj 3 = £34,5 V* V s . On 
vorticity equation, accompanying (|13|) with (JHJ) and taking account of 9 = 0, one has 



LU 



(V x a) • n s = e m V k ai 

}2 



p ax^x 1 p V J p 



Thanks to the relation (see IXie et all . 120131 ) 

e M3 V fe (V s V s ^) = V s V s (e m V fc ^) + e M3 V fe (A G ^) = V s V s u 3 + e m V fc (A G ^) 

the proof is completed. 

Finally, the pressure Possion equation can be derived 

Lemma 7 (Pressure Possion Equation for Incompressible Flows) 

-Ap = p[(V ® V) : (V ® V) + K G | F| 2 ] - 2//V • (VA G ) - V • / s , |F| 2 := F S V S 
Proof: Taking denoted as V- the divergence operator on both sides of (JSJ), one has 

pV l (^(x,t) + V s V s vA = -V l (V l p) + p[v l V s V s V l + V l (K G V l )] +V l f sur<l 

where 

v\v s v a Vi) = (v l v s )(v s Vi) + y s (v'v s ^) = (W s )(v a v*) + y s [v s (v^) + Rf l : s v t 

= (V l V s ){V s V{) + Kg^I - g a g ls )V s V l = (V V) : (V ® V) + ff G |V| : 



V'(V S V S ^) 



VV(V S V*) = V S V Z (V S V^) + i? s " s V t V* + RfV.V t 

V S V Z (V S ^) + A G (^ ts - g tl SI)V t Vi + A G (^ S - g u Sl)V a V t = V s (vV s ^) 

V S [V S (V^) + flftVj] = V S [A G (<^* - g tl g ls )V t ] = V S (K G V S ) = V S V S K G 



Then the identity is proved. 



Fixed Solid Botmchn 





Figure 2: Sketch of local parametrization of a patch of surface, D x i x 2, D x i x 2 and D x » x ± denote 
parametric domains with respect to coordinates {x l ,x 2 }, {x 1 ,^ 2 } and {x'^x- 1 } respectively. 

4 Conclusions 

Basically, some primary relations in general theory of vorticity dynamics have been extended 
to surface tensor fields based on Levi-Civita connection operator. Subsequently, the governing 
equation of vorticity, Lagrange theorem on vorticity, Caswell formula on strain tensor are at- 
tained. Furthermore, the stream function & vorticity algorithm with pressure Possion equation 
has been set up for incompressible flows. It has been revealed that Gaussian and mean cur- 
vatures are taking part directly in some governing equations but all of the effects due to the 
geometry of the surface will disappear as flows on general fixed surfaces degenerates to ones on 
planes. All of the theoretical results presented in this paper are exactly deduced without any 
approximation and constitute a theoretical framework of vorticity dynamics for two dimensional 
flows on fixed smooth surfaces. 
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Appendix 1. Construction of local curvilinear coordinates {x [ \x^} 

As sketched in Figure [51 one can suppose without lost of generality that the preimage of a 
segment of boundary on the surface in the parametric space {x l ,x 2 } can be represented locally 
as 



(a, b) B x 1 H> 



(x 1 ) 



^(x 1 



Then the following coordinates can be constructed locally 

■-1 



x 



x 1 



x 



(x\x 2 ) 



X 



,1 
ijj{x l 



X 
2 J,(„l^ 



The segment of the boundary is corresponding to x 1 G (a, b) and x 2 = 0. The local co-variant 
basis with respect to {x*} 2 =1 is denoted by {<7j} 2 =1 . Thirdly, one introduces another coordinates 
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{x",x }. Its local co-variant basis {g u ,g±} has the following relations with {<7j}f =1 



\9\»9±.\ = [9i,92l 

dx" dx 1 - 
Accompanying the construction of {g n ,g±} on the boundary with the general relations 



dx 1 dx 1 



dx" dx 1 - 
dx 2 dx 2 
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dx 1 dx 2 

dx 1 - dx L 



dx 1 dx 2 



(x\x 2 ) 



7,1 



dx 
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dx" dx 1 - 

dx 2 dx 2 



dx" dx 1 - 



(x",x ± ) 



the following Jacobian matrix on the boundary can be determined 

(s\0), atefab) 



dx" dx" 
dx 1 dx 2 

dx dx 



- dx 1 dx 2 
then the coordinates {x",x^~} can be constructed as follows 



x"(x\x 2 ) 



x (x , X ) 



Sr^.oj^ + ^CxSoj-x 2 



dx 
dx 



dx -(x\Q)dx 1 + % I {x\G) 



dx 2 

dx~ 
dx 2 
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